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We consider the top-bottom doublet in the background of the sphaleron for 
the realistic case of large non-degeneracy of fermion masses, in particular 
rrib = 5 GeV and m t = 175 GeV. We propose an axially symmetric (r, 0)- 
Ph. dependent ansatz for fermion fields and investigate the effects of the non- 

degeneracy on them. The exact solution is described, with an error less than 
0.01%, by a set of ten radial functions. We also propose an approximate 
solution, in the m^jmt —>■ limit, with an error 0(mb/m t ). We have found 
$_i ■ that the effects of non-degeneracy provide a ^-dependence typically ~ 10%. 
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After the discovery by t'Hooft JTJ that baryon and lepton numbers are anomalous 
global symmetries in the Standard Model, violated by non-perturbative effects, and the 
observation |2j that the three Sakharov's conditions for baryogenesis can be fulfilled at 
the electroweak phase transition, there has been renewed interest in understanding the 
dynamics of baryon asymmetry generation near the electroweak critical temperature. 
In particular the sphaleron, a static unstable solution of the classical equations of 
motion H, has played a key role because the baryon violation rate for fluctuations 
between neighboring minima contains a Boltzmann suppression factor proportional to 
the energy of the sphaleron £7 sp h- In particular the rate per unit time and unit volume 
is given by || 

T~ujT exp I - — 1 (1) 

where the sphaleron solution is usually computed in the limit of zero Weinberg angle 
(corrections of O(g') have been proved to be negligible || in the absence of fermions) 
and use has been made that the energy of the sphaleron at finite temperature follows 
an approximate scaling law with (f)(T)/v |J. The prefactor u in (|IJ) which contains, 
in particular, the product of bosonic and fermionic fluctuation determinants (ftbos^fer) 
around the sphaleron, was evaluated J7| in the high T limit (three-dimensional effective 
theory with all fermionic modes decoupled, i.e. Kf er = 1). 

In view of the recent experimental detection of the top-quark and its mass mea- 
surement H it seems of utmost importance to consider the modification, due to the 
presence of the top-quark, of the classical sphaleron energy and its contribution to /tf er . 
Progress in this direction has been done by considering both effects, in Refs. and 
T0[| , respectively, as triggered by the presence of mass degenerate fermion doublets. 



The reason for studying degenerate fermions is simplicity since the spherical symmetry 
of the sphaleron solution is not spoiled by the presence of degenerate fermion doublets. 

However, for the top-bottom doublet, with hierarchically different masses, rrif, ~ 5 
GeV and m t ~ 175 GeV, the breaking of degeneracy (custodial symmetry), controlled 
by (m t — mi,) / m t , cannot be considered a priori as a small perturbation and the previous 
results might be spoiled by the large non- degeneracy effects. It is the purpose of 
this note to study the effects of the non-degeneracy of the top-bottom doublet in the 
background of the sphaleron. Our results have to be considered as a first step towards 
quantifying the effects of the non-degeneracy of fermion doublets on the sphaleron 
energy. They contain interesting (and somewhat unexpected) results which could be 
used as hints towards computing the effects of the top-bottom doublet with physical 
masses on the classical sphaleron barrier and fermion determinant. 



After completion of this work we have received Ref. UTTJ where the same problem 
is considered using a different expansion for the fermion fields. We will compare both 
approaches and point out the differences in both expansions. 

We will now study the fermionic distribution of a realistic quark doublet, non- 
degenerate in mass, in the sphaleron background. To simplify, we will work in the 
g' = approximation. In this limit the hypercharge field decouples and the 
sphaleron energy density is spherically symmetric. The sphaleron configuration can 
be parametrized as: 
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W? = f r e ai] r 3 f(r) 



The relevant lagrangian density for fermions is: 

C = iQ$Q + iu R $u R + id R $d R 



where 



-h t (Q®u r + h.c.) - h b (Q$d R + h.c. 

Q = 



(3) 




is the left-handed quark doublet, = d^W^-d^W^+ge^W^W^ = d^-i^W", 
The equations of motion for the fermionic fields can be written as: 



z^Q = h t $u R + h b <$>d R 



(4) 



The non degeneracy of fermions explicitly breaks the global custodial symmetry 577(2) l x 
SU(2) R . Let us remember that in the degenerate case the Dirac hamiltonian deduced 
from (|4]) commutes with the grand spin operator K = L + S + T, where L stands for the 
angular momentum, S for the spin and T for the isospin. In this case the eigenstates \I/ 
of the Dirac operator can be labelled by k, k s , with K 2 ^> = k(k + and K^> = k^. 
In particular the zero mode is included in the k = subspace. For non-degenerate 
fermions, the different Yukawa couplings to the Higgs field background single out a 
direction in the isospin space. As a consequence K 2 does not commute any longer with 
the hamiltonian and the energy eigenstates will have non-vanishing projections over 
different values of k. However, K 3 still commutes with the hamiltonian, and hence k 3 
is a good quantum number to label the energy eigenstates. Since k 3 is a discrete label 
and we want to continuously deform the zero mode from the degenerate case to the 
non-degenerate one, we will restrict ourselves to the subspace k 3 = 0. 
To describe the fermions, we deal with the Dirac representation: 
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with a 1 the Pauli matrices and a = 1. Using two component spinors, qi, q R , defined 
as: 
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we can write our ansatz, a general parametrization of the k 3 = subspace, as: 



u L (r) 


= (u L2 (p,z) +u L1 (p,z) p.o)x- 




— \ULl\P, z) -r UL2\P, Z) p. a) x+ 


u R {r) 


= (u R2 (p, z) + u R1 (p, z) p.a) x- 


d R (r) 


= (d R1 (p, z) + d R2 (p, z) p.a) x+ 



(7) 



where we have decomposed r into the longitudinal and transverse components, f = 
p + zk, p is the normalized vector in the XY plane p/p, p.a = p x a x + p y a y and 
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The equations for fermions can be written as functions of r, 9, with p = r sin 9, as: 



d r u L1 (r, 9) + ±u L1 (r, 9) + ±d e u L2 (r, 9) 

+i/(r) (| sin 29(d L1 (r, 9) + u L2 (r, 9)) + sin 2 9(d L2 (r, 9) - u L1 (r, 9))) 

-h(r) (± sin 29(h b d R1 (r, 9) + h t u R2 (r, 9)) + cos 2 9h t u R1 (r, 9) + sin 2 9h b d R2 (r, 9)) = 

d r u L2 (r, 9) - cot 9\u L1 (r, 9) - }d e u Ll (r, 9) 

- l r f(r) ((d L1 (r, 9) - u L2 (r, 9)) + cos 2 9(d L1 (r, 9) + u L2 (r, 9)) 

+±sin20(d L2 (r,0)-u L1 (r,0))) 
-h(r) (| sin 29(h t u R1 (r, 9) - h b d R2 (r, 9)) + sin 2 9h b d R1 (r, 9) - cos 2 9h t u R2 (r, 9)) = 

d r d L1 (r, 9) + cot 9\d L2 [r, 9) + \d e d L2 [r, 9) 

- l -f{r) ((u L2 (r, 9) - d L1 (r, 9)) + cos 2 9(d L1 (r, 9) + u L2 (r, 9)) 

+ ±sm29(d L2 (r,9)-u L1 (r,9))) 
-h(r) (| sin 29(htU R1 (r, 9) - h b d R2 (r, 9)) + sin 2 9h t u R2 (r, 9) - cos 2 8h b d R1 (r, 9)) = 

d r d L2 (r, 9) + ld L2 (r, 9) - ±d e d L1 {r, 9) 

-i/(r) (| sin 29(d L1 (r, 9) + « L2 (r, 9)) + sin 2 £(d L2 (r, 0) - u L1 (r, 9))) 

+h(r) (\sin28(h t u R2 {r,6) + h b d R1 (r,9)) - sin 2 9 h t u R1 (r, 9) - cos 2 9h b d R2 (r, 9)) = 

d r u R1 (r, 9) + I Ufll (r, 0) + i<9 eMi?2 (r, 0) 

-/i(r)/i t (| sin 20(d L i(r, 0) + w L2 (r, 9)) + cos 2 0u L1 (r, 0) + sin 2 9d L2 (r, 9)) = 
d r u R2 (r, 9) - cot 9\u Rl (r, 9) - \d e u R1 (r, 9) 

-h{r)h t (| sin 29(u L1 (r, 9) - d L2 (r, 9)) + sin 2 9d L1 (r, 9) - cos 2 ^ L2 (r, 9)) = 
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d r d m (r, 9) + cot 9±d m (r, 9) + \d e d R2 {r, 9) 

-h(r)h (§ sin 29(u L1 (r, 9) - d L2 (r, 9)) + sin 2 9u L2 (r, 

d r d R2 (r, 9) + id R2 (r, 9) - ±d e d m (r, 9) 
+h(r)h b ^sm29(u L2 (r,9) + d L1 (r,0)) - sin 2 9u LX (r, 



— cos 2 Odnir, 



cos 2 9d L2 (r, 



(9) 



We have, then, a set of eight partial differential equations. Instead of solving it in 
general, we will work out a systematic expansion of these functions in terms of a Fourier 
series. Notice that the commutator of the hamiltonian [defined on the (q^, g^) T -space] 
with parity vanishes on the zero modes, which allows to label them with a definite 
parity. As parity is a discrete symmetry we will choose for our ansatz the same parity 
(i.e. even) as zero modes in the degenerate case. Therefore, our expansion can be 
written as 
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We expect the convergence of the series ( [Tof) to be good enough to provide a reliable 
solution to the system of equations (Jffy. In fact, cutting off the series at n = would 
include the appropriate ansatz to describe a situation in which the fermion doublet is 
degenerate. However the case we want to address, i.e. that of the top-bottom doublet, 
is very massive and highly non-degenerate. This means, on the one hand, that we need 
the modes n > 1 in ([TDD an d ; on the other hand, that the convergence of the series 
(|10|) might be slow enough to require many terms for an accurate description of the 
solution. Fortunately we will see this is not the case, and, as we will prove, taking 
only the terms with n — 0, 1 provides a reliable solution with an accuracy better than 
0.01%. 



1 An alternative expansion is given by developing the fermionic fields in terms of the even spherical 
harmonic subsets Y2 n (6), Y£^(9,<f>) instead of cos(2rt6>), sin(2n0)e ±I< ^ respectively. In this basis it is 
straightforward to check that = 0. 
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Assuming the convergence of series fllUD one can solve the system (^j) to any given 
order n (i.e. neglecting all modes with higher values of n) and estimate the accuracy 
of the different approximations. A first step along this direction is solving (|9]) to order 
n = 0, i.e. in the approximation of a spherically symmetric ansatz for fermions (as in 
the case m u = m^). Only four spinor components do not vanish: 

u (°) d {0) u {0) d {0) (11) 

a L2i a Lli U R2' U R1- \ LL J 

Replacing them into (Q) one obtains a system of four differential equations whose 
solution is shown in Fig. 1 (dashed lines) 0. It is interesting to notice that the relation 

«S(0 + dS(r) = 0, (12) 

valid in the degenerate case, also holds in this approximation. On the other hand, 
the similar relation valid in the degenerate case, i4m r ) + ^m( r ) = is spoiled by 
non-degeneration effects, and replaced by 

h t d%l(r) + h b u%l{r) = 0, (13) 

In fact, the functions defined by the left-hand sides of Eqs. (|l^) and (|13|) only couple 
to themselves and to higher modes, and hence can be consistently fixed to zero to this 
order. 

In general, solving the system to a given order n, amounts to keeping the system 
of fields 

7 >) 7/ (p) M M 

a Lli a L2i u Lli u L2i 



7 >) 7 >) JPi f]<P) (1A) 
a Rl: a R2i a R\i a R2 \ Lq: ) 

for p — 1, ...,n, and the spherically symmetric components fllTI): a system of 8n + 4 
differential equations with 8n + 4 unknown functions. We will solve the system of 
differential equations using the following boundary conditions. There are two zeroes 
in the (8n + 4) x (8n + 4) asymptotic mass matrix which correspond to two directions, 
in the space of highest, n-th, components of the spinorial fields. They must be fixed 
to zero at infinity, and thus elsewhere since these directions have (up to higher order 
contributions) fixed points at zero. Fixing to zero these two directions amounts to the 
conditions: 

u£?(r) + u&>(r) + dg(r) - cfeV) = (15) 

and 

MS(r) + Mg(r) + Mg(r) - Mg(r) = (16) 

Out of the remaining 8n + 2 fields, An + 1 correspond to positive mass eigenvalues, 
and therefore are fixed to zero at infinity imposing normalizability of the fermionic 
configuration, and 4n + 1 have negative mass eigenvalues. An of them are fixed to zero 



2 All dimensional quantities are expressed in units of the corresponding power of v = (<^o)) and we 
have arbitrarily normalized the total fermionic density to 47r. 
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at the origin (triggering the vanishing at the origin of all 8n ^-dependent components), 
and the last one is used to fix the global normalization. In particular, for n = 1 we 
obtain a system of twelve differential equations with twelve unknown functions: the 
set (PD and ([14]) with p — 1. The solution is presented in Figs. 1 and 2 (solid lines). 
We can see from Fig. 1 that the relations (^) and (13) are spoiled by the presence in 
the equations of the angular dependent components. The spoiling is ~10%, which is 
the expected order of magnitude of the angular dependent part of the solution. The 
latter, i.e. the radial functions of ( |14"D with p = 1 are plotted in Fig. 2 (solid lines). 
We can see that the typical size of these functions is ~ 10% those contributing to 
the spherically symmetric part of Fig. 1. This is the amount by which the relations 
between the components with n = 0, Eqs. fll~2|) and (|13"D , are spoiled. The latter are 
replaced by the similar relations ( |I5[ ) and ([16]) between the n = 1 components. One can 
easily check from the curves in Fig. 2 that relations (]15|) and (|16|) are exactly satisfied. 
Moreover, we have numerically found that Eq. ([15]) splits with a good accuracy, for 
n > 1, into the couple of constraints: 



r) + d 
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(18) 



We believe that Eqs. (|17D and (|1^) are not a numerical coincidence. The fermionic 
fields expanded up to a given order n can be decomposed into eigenstates of angular 



momentum 



0, 2, 



< 2n, isospin T3 = ±1/2 and spin S3 = ±1/2, given by 



\£,-(T 3 + S 3 )) 



Defining J = L+T and using the basis of J 2 , J3 eigenstates, condition ( |l7f ) is equivalent 
to projecting out the left-handed states on \j, j 3 ) = 2n± |, — |V Similarly, condition 

2n — I , h > ■ Therefore a consistent ansatz 



(pTq) amounts to projecting out on \j, j s ) 
using the latter basis should take advantage of this approximate symmetry. 

We have solved the system (|9]) to order n = 2. In this case we have a system of 
twenty differential equations with twenty unknown functions, those in (111) with p = 1,2 
and in (pTTj) . The new components, with p = 2 in ([14]), are less than 10~ 4 the large 
components in Fig. 1 while the p = 1 components in ( |T4| ) and the components in (|TT[), 
computed to this order of approximation, are indistinguishable from those previously 
computed. In fact, we have plotted in Figs. 1 and 2 the latter in solid lines and they 
are superimposed on those computed to order n = 1. In the same way the relationships 
(|17|) and (|18[ ) remain very accurately satisfied. 
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To summarize the previous results, we can keep, out of the general expansion fllCf), 
the particular reduction 



u L (r) 



u Li(. r ) sin2#e 



U L2 ( r ) cos 2^ 



u R {f) 



u 



(0), 
ii'2 ! 



"4i( r ) sin2#e 



r) + m^2(^)cos2^ 



41(0 



+ d l H(r) cos 20 



4 1 2 ) (r)sin2^ 



41(0 



+ 



dg(r)cos20 



d^(r)sin2^ 



9) 

where two radial functions can be eliminated using the relations fll5|) and (|l6l). 

Given the smallness of the ratio m^/nit one can try to implement the mb/rrit —>■ 
limit in our solution. In that limit d R decouples from the equations of motion and can 
be consistently fixed to zero. In particular 



41(0 



41(0 



= 0. 



(20) 



Eq. ( P0|) is very accurately satisfied, as can be checked from Figs, lb and 2b. We have 
solved the problem for values of m& in the range < < 5 GeV. We have verified in 
this range that d R is proportional to h b - The rest of fields remain essentially unchanged 
since their dependence on mj, is through 0(til) terms. Numerically these non- vanishing 
fields, computed in the strict limit m b = 0, are superimposed with the corresponding 
functions plotted in Figs. 1 and 2 for the physical value of nib- We can conclude that 
neglecting all right-handed components of spinor fields, in the limit vn^jmt — > 0, is 
accurate to 0(m b /m t ). 



We will compare now the previous solution with that recently proposed in Ref. [1 1 
The solution in Ref. [[□]] corresponds to the parametrization in terms of eight radial 
functions, as: 



u L (r) 



- \F+ sin 26e 



Gt + \Ft+ 



|F/cos2# 



u R (r 



dr(f) 



±F+ sin 29e 



\F++ ±F+cos2fl 



\FZ+ \F£cos29 



\F^ sin 29e^ 



(21) 



d R {r) 



-G R + \F R -+ \F R cos26 



\F R sin26e^ 
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Ff = FAr) ± AFAr), G$ = GAr) ± AFAr) 



where 

(22) 

F R = F R (r) ± AF R (r), G± = G R (r) ± AF R (r) 



Since we have proved that our expansion ([19]) describes the exact solution with 
an error less than 0.01%, we can compare it with the parametrization in Eq. fl2ip . 
In the left-handed sector Eq. ( f2~T|) would imply that «£i + u^l — and $l\ = d^: 
these conditions are consistent with ( |T3| ) though a quick glance at Fig. 2 shows that 
they are violated in the exact solution by ~10% of the corresponding components, an 
error much greater than that induced by the next order. Moreover, these conditions 
conditions are not controlled by m^/rrit but, contrarily, become exact in the limit 
of large degeneracy. Similarly, in the right-handed sector, Eq. ([H]) would predict 
u^ R1 + u R l = and d R \ = d R2 , which are also consistent with our n — 1 constraint 
(|16|) but not satisfied by the exact solution. In summary, we have found that the 
parametrization (21) is accurate in the case of large degeneracy, but has corrections not 
controlled by the parameter mb/m t , which are much greater than those corresponding 
to the next order in our expansion. 

In conclusion, we have performed a consistent expansion for the top-bottom fermionic 
fields in the background of the sphaleron. We have found that the first two terms in the 
expansion provide the fermionic fields with an accuracy better than 0.01%. We have 
quantified the effect of the top-bottom non-degeneracy in the ^-dependent part of the 
fermionic components as ~10%. We have proved that the approximation m^jnii — * 
provides a solution accurate to 0(rrib/m t ). Our work should be considered as a first 
step towards evaluating the effects of the non-degeneracy of the top-bottom doublet in 
the sphaleron energy. 



8 



References 



G. t'Hooft, Phys. Rev. Lett. 37 (1976) 8; Phys. Rev. D14 (1976) 3432 

V.A. Kuzmin, V.A. Rubakov and M.S. Shaposhnikov, Phys. Lett. B155 (1985) 
36; Phys. Lett. B191 (1987) 171 

N.S. Manton, Phys. Rev. D28 (1983) 2019; F.R. Klinkhamer and N.S. Manton, 
Phys. Rev. D30 (1984) 2212 

P. Arnold and L. McLerran, Phys. Rev. D36 (1987) 581; Phys. Rev. D37 (1988) 
1020; S.Yu Khlebnikov and M.E. Shaposhnikov, Nucl. Phys. B308 (1988) 885 

J. Kunz, B. Kleihaus and Y. Brihaye, Phys. Rev. D46 (1992) 3587 

Y. Brihaye and J. Kunz, Phys. Rev. D48 (1993) 3884 

L. Carson, X. Li, L. McLerran and R.-T. Wang, Phys. Rev. D42 (1990) 2127; 
J. Baacke and S. Junker, Phys. Rev. D49 (1994) 2055 

F. Abe at al, CDF Collaboration, Phys. Rev. D50 (1994) 2966; Phys. Rev. Lett. 
73 (1994) 225 and preprint FERMILAB-PUB-95/022-E; 

S. Abachi et al., DO Collaboration, Phys. Rev. Lett. 74 (1995) 2422 and preprint 
FERMILAB-PUB-95 /028-E 

J. Kunz and Y. Brihaye, Phys. Lett. B304 (1993) 141; G. Nolte and J. Kunz, 
Phys. Rev. D48 (1993) 5905 

D. Diakonov, M. Polyakov, P. Pobylitsa, P. Sieber, J. Schaldach and K. Goeke, 
Phys. Lett. B336 (1994) 457; Phys. Rev. D49 (1994) 6864; and, St. Petersburg 
preprint |hep-ph/9502245| 



G. Nolte, J. Kunz and B. Kleihaus, Oldenburg preprint [hep-ph/ 9507369 



9 



0.75 



0.50 - >! 



0.25 - 



0.00 



0.75 



0.50 



0.25 



0.00 




Figure 1: a) The radial components U L 2 and Uj^ computed to order n = (dashed lines) and 
n = 1 and 2 (solid lines), b) The radial components — d^i and —d^\ computed to order n = 
(dashed lines) and n — 1 and 2 (solid lines). 
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